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Radiative El decay widths of X(3872) are calculated through the relativistic Salpeter method, 
with the assumption that X(3872) is the Xci(2P) state, which is the radial excited state of Xci(lP). 
We first calculated the El decay width of Xci(lP). The result is in agreement with experimental 
data excellently. Then we calculated the width of X(3872) with the assignment Xci(2P). Results are: 
T(X(3872) -> 7J/VO = 33.0 keV, P(X(3872) ->■ 7^(25')) = 146 keV and T(X(3872) -> 7^(3770)) = 
7.09 keV. The ratio Br(X(3872) ->■ 7i/>(2S))/Br(X(3872) -» 7J/VO = 4.4 agrees with experimental 
data by BaBar, but is larger than the new up-bound reported by Belle recently. With the same 
method, we also predicted the decay widths, r(xbi(lP)) -> 7T(1S)) = 30.0 keV, r(x6i(2P)) -¥ 
7T(1S)) = 5.65 keV and r(xw(2P)) -> 7T(2S)) = 15.8 keV, from which we get the full widths: 
r(xt>i(lP)) ~ 85.7 keV and r( X 6i(2P)) - 66.5 keV. 

PACS numbers: 13.25.Gv, 13.40.Hq, 13.60.-r, 13.25.Ft 



I. INTRODUCTION 

X(3872) was discovered by Belle Collaboration |l| in 2003 through the channel B ± -> K ± J/^tt+tT '. The mass 
reported by Belle is M = 3872.0 ± 0.6 ± 0.5 MeV, and the full width has an upper limit T < 2.3 MeV at 90% C.L.. 
Later the existence of this particle was confirmed by CDF [2], DO [jj, and BaBar 4| Collaborations. 

The radiative decay channel X(3872) — > J /ipl |5J indicates that this particle has positive C-parity. The most 
possible J PC of X(3872) is 1 ++ , which is favored by the analysis of the decay angular distribution [6j. But the dipion 
mass distribution and tripion mass distribution in X(3872) — > J/ipn + Tr~ and X(3872) — ► J / ipnr + n~ ir° [5| favor a p 
resonance and a ui resonance, respectively, which indicates a large isospin breaking. The mass of X(3872) is 50-100 
MeV smaller than the predictions of potential models. To understand these puzzles, many assignments of X(3872) 
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were propose d, besides the traditional charmonium state assignment 7 
state 



lfj , there are the assignments of a molecular 
18j ]. cusp effect or virtual state Q, [3| (for a 



16| . a hybrid charmonium [17J, a diquark-antidiquark state 
review, see e.g. Ref. |21j). 

If X(3872) is a 1 ++ state, then there are two most possible natural assignments, a molecular state or a traditional 
charmoniumx c i(2P). Molecular state model predicts the value of Br(B+ -> X(3872)K+)/Br(B° ->■ X(3872)K°) is 



about 10% 



22} , while the experimental value is 0.5 ± 0.3 ± 0.05 



used to calculate radiative decays. From Ref. 
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231 ] . This model meets more serious problems when 
(see Table. 1), one can get the ratio r^,' 7 /rj/^, 7 ~4x 10~ 3 , while 



the experimental value by BaBar is 3.4 ± 1.4 24]. In this Letter, we will not consider the possibility of X(3872) as a 
molecular state, but due to its El radiative decay, we consider the possibility of an ordinary charmonium state. 
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Because El radiative decays will play a fundamental role in determination of the nature of X(3872), in this Letter, 
we just calculate the radiative El decay widths of X(3872) by assigning it as the Xci(2P) state and give the results. 
Although there is a discrepancy in the mass values of experiments and models, as Ref. [8j proposed, this is due to 
additional effects, such as coupled-channel effect. For the large isospin breaking, charmonium model can also give a 
good explanation 25|. 



This Letter is organized as follows. In Sec. II, we solve the instantaneous Bathe-Salpeter (BS) equation (Salpeter 



tion) flfl, 



equation) [26l l27j . and get wave functions of initial and final states. Then within Mandelstam formalism 
calculate the transition matrix element. In Sec. Ill, we compare our results with other theoretical predictions and 
experimental data, some predictions and discussions are also given in this section. 

II. El DECAY OF X(3872) WITH Xci(2P) CHARMONIUM ASSIGNMENT 

The wave function of 1 ++ sate is , 

*»i+(gx) = ie^ aP p v qi4[viMr + <p 2 pr + + viPr4JM\/M 2 , (i) 

where e^ a p is the totally antisymmetric tensor. t\ is the polarization vector of the meson while M is its mass. P and 
q are the total momentum and relative momentum of constitute quark and antiquark, respectively, which are defined 
as: 

mi TO2 

Pi = a x P + q, ati = ■ , p 2 = a 2 P — q, a 2 = — , (2) 

mi + m 2 mi + m 2 

where p\,p 2 are the momenta of quark and antiquark, respectively, mi = m 2 is the mass of constitute quarks, (pjS 
are functions of q\. q±_ has the form: q 1 ^ = q 11 — {P ■ q/M 2 )P^ . Because there are two constrain conditions 29j, 
tp3, (^4 can be expressed by <pi, tp 2 [29]. The wave function above has a different form with that in 29], but they 
are equivalent to each other. We show a general wave function form for 1 + state, which means quark and antiquark 
inside the meson can have different masses. If we consider charmonium 1 ++ state, the quark and antiquark have the 



same mass, then <p 3 will disappear 29]. 
The wave function of 1 state is 3fJ, 

H>i-Wx) = l'±^2 [fi(q'±) + P f f 2 (q' ± )/M f + jM ± )/M f + Pfff x Mq' x )/M*] + M f &Mq' ± ) 

+ &Pfh(q'x) + iAi 2 - q'x • *2)frtf±) + {Ptbh - Prfx ■ e 2 )h(q' x )/M f , (3) 

where Mf, Pf and e 2 are the mass, momentum and polarization vector of the meson, respectively. Again, if we 
consider charmonium, the constitute quark and antiquark inside the mason have the equal mass. Because there are 
four constrain equations [30(, /V and f 2 will disappear, and fx and fs can be expressed by f^, f$, f$ and fe- Here we 



will not present the details of solving BS equation, which can be found in Ref. 



3l|. 
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FIG. 1: Feynman diagram for the transition: \ c i — > J/^/j + 7. 

We just give the Cornell potential which is applied when solving BS equation: 

V(q) = V s (q) + l0 ® 1 °V v (q), 



V s (q) = -(- + V )5 3 (q)+ X ' 



Q 



r2 ^2 



(<7 +a 2 y 



V v (q) = - 



2 a„(g) 
3 ^ 2 7f ^ 2' 



a s (9) 



12tt 



25 



ln(a 



Aqcd 



(4) 
(5) 

(6) 
(7) 



Here A, a, e, Vb and Aqcd are parameters. By fitting the mass spectra of 1 ++ , 1 masons, we can find the best-fit 
values of these parameters: a = e = 2.7183, a = 0.06 GeV, A = 0.2 GeV, m c — 1.7553 GeV, m& = 5.13 GeV, 
Aqcd = 0.26 GeV (cc), 0.20 GeV (66) (see Q). For 1++ state, V =-0.452 GeV (cc), -0.521 GeV (66), for 1~ 
state, Vq — -0.465 GeV (cc), -0.570 GeV (66). Here a is the effective gluon mass. Since the potential we chose is a 
phenomenological one and the gluon mass as a parameter is not running, the value of a here is lower than the usual 
chosen especially when it is running close to the infrared limit. 

Wave functions above are constructed based on the quantum number J p or J PC of mesons. For example, J p of 
every term in Eq. (3) is l - (or 1 for equal mass system). One can see that there is S wave and D wave mixing 
automatically, especially for the third state (?/>(3770)), which is D wave dominating, but mixing with a small part of 



S wave. This can be seen clearly in spherical polar coordinates 

The relativistic transition amplitude of 1 ++ state decaying to a photon and a 1 state (see Fig.l) can be written 
in terms of BS wave function: 

{2Tr) 4 ee q 



T= (P /e2 ,fce|S|P ei ) = 



-S 4 (Pf +k- P)e 6 M £ , 



(8) 



where e, t\ and e 2 are the polarization vectors of the photon, initial meson and final meson, respectively. P, Pf and 
k are the momenta of initial meson, final meson and photon, respectively. e q = | for charm quark and e q = — i for 
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FIG. 2: Radial wave function of X(3872) and f5 of J/tp. 




FIG. 3: Radial wave function of X(3872) and f5 of ^(2S). 
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FIG. 4: Radial wave function ^-<pi of X(3872) and J^-f3 of -0(3770). 
bottom quar k are the charges in unit of e. is the matrix element of the electromagnetic current, which according 



to Refs. 



28 



32j, in the leading order (the order of a — f-, also neglect terms contain ip + , ip + and ip > which 



contribute less than 1%) can be written as: 

— > 

where tp ++ is the positive part of BS equation. ftj_ and (p ++ are defined as Pj ± = Pj — (P ■ Pf /M 2 )P fl and 
7o(</? ++ ) t 7o, respectively. 

For X(3872), the positive energy part of wave function has the form: 

= ie^P v q a A{A x r + A a7 V + 4j7V*l)> (10) 

where A\, A 2 , A% are defined as: 

^ = 2^ + mB ) «' (12) 
2M m M ui 

The positive energy part of wave function for 1 state can be written as: 

p+t = Si ^2 + S 2 £ ft- + B 3 ft / 2 / ± + S 4 g ± • e 2 + S 5 <zi ■ e 2 ft + S 6 <?1 ■ e 2 ^ + S 7<Z± • e 2 ft- (14) 



where the expressions of B\ to B-j can be found in Ref. 



TABLE I: El decay widths of Xci(lP) an d Xci(2P). In Ref. Q, Barnes and Godfrey (labeled by B&G) have made the impulse, 
nonrelativistic, zero recoil, and dipole approximations. In Ref. [T3|, three method are adopted: one has the same approximations 
as B&G, but compute with a improved potential (labeled by Swansonl); one has no approximation (labeled by Swanson2); the 
last one is molecular model (labeled by Swanson3). Our values inside the parentheses are for the cases that the mass of 3923 
MeV for X(3872) is chosen. 



1 1 r 

Her. 


r j/V 7 


r W 7 ( keV ) 


r V-(2S) 7 ( keV ) 


^(STTOVr ( keV ) 


This work 


306 


33.0 (33.3) 


146 (182) 


7.09 (9.83) 


PDG [34] 


320 








Li and Chao [7] 




45 


60 




Swansonl [11] 




71 


95 


6.5 


Swanson2 [11] 




139 


94 


6.4 


Swanson3 [11] 




8 


0.03 





B&G [8] 




11.0 


63.9 


3.7 


Eitchen et al [35] 




110 


180 


25 


Dong et al [36] 




1 ~ 2 


5 ~ 6 





Ref. 



TABLE II: El decay widths of Xm(1P) and x&i(2P). 



jfgg (keV) 



rffi s ff (keV) 



Y(2S) 7 



(keV) 



This work 
Kwong and Rosner [38J 
Ebert et al [39] 
Fazio [9] 



30.0 
32.8 
36.6 
107 



5.65 
9.31 
7.49 



15.8 
15.9 
14.7 



III. NUMERICAL RESULTS AND DISCUSSIONS 

We first calculated the decay width of Xei(LP) to J/i/j and 7. The result 306 keV shown in Table. 1 agrees with 
the experimental value 320 keV very well. This shows that our method can be used to describe radiative decay. For 
X(3872), with the 2 3 P X charmonium assumption, we calculated decay widths of three channels. We first solved the 



29(, which 



instantaneous BS equation by setting the parameter Vq — —0.452 GeV. The mass of Xci(2P) is 3.923 GeV 
is about 50 MeV larger than that of X(3872). This is the common character of all potential models, which may be 
due to the coupled-channel effect. The results which we got by using this wave function are included in parentheses 
in Table 1. To make the mass of Xci(2P) equal to 3872 MeV, we solved the BS equation by setting Vb = —0.516 GeV, 
and keeping other parameters un-changed. (This also causes a mass decrease of 50 GeV for other states. Here we just 
want to get the wave function of x c i(2P) when its mass is 3872 MeV. To make the spectrum agree with experimental 
data, we have to modify our coupled equations, especially the potential, which is our future work.) Decay widths 



with this set of parameters are those outside parentheses. We can see that the value of r is almost unchanged, 

tively. 



J/V>7 

while the values of F^icfv? an< ^ ^iT^rrm-w are reduced by nearly 20% and 30%, respec 



One can see that our V = 33.0 keV is of the same order with that of Li and Chao 



*1 



B&G 



and Swansonl 



21 



but much larger than 8 keV of Swanson3 1.1] (molecular model) and 1 ~ 2 keV of Dong 1 36] (molecular and cc mixture) 
Our r^/ 7 = 146 keV is about 2.5 times larger than that of Li and Chao 7] and B&G [8j, but approximately equals to 
that of Eitchen 35] , which has considered the influence of open-charm channels. The results in Swanson2 [111 ] have used 
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and B&G 



a improved potential and included no zero recoil and dipole approximation which used in Swanson3 
But as B&G Q did, the wave function and meson mass are calculated by adding spin-dependent interaction in the 
Hamiltonian. In this Letter, we started from BS equation, which is relativistic covariance. By using instantaneous 
approximation, we get coupled Salpeter equations, which has included the relativistic effects automatically. 
The ratios of El decay widths and the width of X(3872) — > 7r + 7r~ J/tjj detected by BaBar are 

Br(X(3872) -»• 7 J/V0 



2J: 



Br(X(3872) ->• n+ir-J/ip) 
Br(X(3872) -> ~fip(2S)) 



= 0.33 ±0.12, (15) 



= 1.1 ±0.4. (16) 



Br(X(3872) -> n+TT-J/ip) 

Up to now, the widths and branch ratios of El decay channels have not been measured precisely. But the ratio can 
be drawn from Eas.lfTB]) and Eq.® |24|: 

Br(X(3872) — > 7?/>(2S)) _ 

Br(X(3872) -> 7J/VO V ' 

With our results in Table 1 we get this ratio: 

Br(X(3872) ^(2S)) 

Br(X(3872) -> -yj/tp) ' v ' 

n n 

which is very close to that of Eq. p7|) . In Refs. [7j and [8j this ratio is 1.3 and 6.1, respectively. We can see that 
models with charmonium assumption can predict this ratio correctly, while molecular model prediction is very small. 
In Ref. 



a composite state which contains both molecular hadronic component and a cc component was considered. 
By changing the mixing angle, a correct ratio can be reached, but the decay widths are dramatically changed. 

Recently Bhardwaj reported the new results of Belle on X(3872) at the QWG2010 conference [37J, which is 
r (x ^3872 )"-»tJ/ ^ ) ^ < 2.1. Our result with the % c i(2P) assignment is two times larger than this up-bound, so there is 
still long way to go to know the nature of X(3872). 

The large ratio r^< 7 /rj/^ 7 can be understood by Figs. 2 and 3. For J/ip, its wave function has no node, that is the 
numerical values of the wave function in the whole space are all positive (see Fig. 2), while for X(3872) and ip(2S), since 
they are the radial excited states of Xci(lP) and J/ip, respectively, the wave functions have one node, that is, before 
the node the values of wave functions are positive, after the node the values are negative. So when we calculate the 
transition amplitude, we need to compute the overlap integral shown in Eq. @. There exists dramatically cancellation 



8 

in the overlap integral before the node and after the node when we consider the decay X(3872) — > jJ/tp which can be 
seen from Fig. 2. This is the reason why the decay width 33.0 keV of this channel is much smaller (almost one order) 
than the width 306 keV of channel Xci(lP) — > jJ/tp. But for the decay X(3872) — > jip(2S), the two overlapping wave 
functions both have the node structures, see Fig. 3. So only in the region where one wave function is before the node, 
while the other is after the node, the overlapping integral gives negative contributions. And we can see from Fig. 3 
that only a very small part of phase space will give negative contributions, so there is almost no cancellation when we 
calculate the transition amplitude. Finally we get a large decay width 146 keV for the channel of X(3872) — > -ftp (2S). 

We have mentioned that the numerical values of r^S)-^ and r^g^gj are very sensitive to the mass of X(3872) 
(see Table 1). This can be explained by different phase space and the node structure of wave functions. From Eq. © 
we can see that in the overlap integral the relative momentum <Zj_ of final state has a shift a.i Pf or —a\ Pf. When 
we change the mass of Xci(2P) from 3872 to 3923 MeV, the node position in the wave functions has almost no change, 
but the value of | Pf | will change obviously due to different phase space, for example, from 181 MeV to 230 MeV 
f° r Xci(2P) — > ^(28)7, which also changes the overlap integral. Finally we got much different values of decay width 
r^^s 2 )^- Similar conclusion can be obtained for the case of Xei(2P) — > 7^(3770) (see Fig. 4). But for Xci(2P) — ► jJ/ipi 
the relative small mass of J/4> results in similar values of | Pf | for both cases, 695 MeV and 736 MeV, so the decay 
widths are similar for both cases. The two-body decay width can be written as: T — ^^j^p-T,\T\ 2 . So the pure 
change caused by the change of phase space is 23.8% for Xci(2P) — > ^/>(2S)7 and 3.2% for Xci(2P) — > yJ/ip. From 
Table 1 the total change of the two processes is 24.7% and 0.9% respectively, which means most of the change for 
Xci(2P) —> "0(2S)7 comes from phase space while for x c i(2P) — > "fJ/ip the larger contribution comes from the change 
of matrix element. 

Using the same method, we also calculated the radiative El decay widths of xm(1P) and x&i(2P), and we show 
the results predicted by our method and other models in Table. 2. One can see that the decay width r(x&i(lP) — > 
7T(1S)) = 30.0 keV calculated by our method is about 3 times smaller than that of Refs. 9J, but close to the values 
in Rcf. 



38| and Ref. 



|, which are 32.8 keV and 36.6 keV, respectively. There are still no experimental data of 
these radiative decay widths, however, ratios are available. Particle Data Group [34[ has listed the branching ratios: 
Br(xbi(lP) -»• 7T(1S)) = (35 ± 8) x 10" 2 , Br(xbi(2P) -> 7 T(1S)) = (8.5 ± 1.3) X lO" 2 , Br(xbi(2P) -> 7 T(2S)) = 
(21 ± 4) x 10~ 2 , so from this experimental data, we can get the ratio [9J: 

Br(xw(2P)^ 7 T(2S)) 



Br(x W (2P)^ 7 T(lS)) 
Our result is 

Br(x W (2P)^ 7 T(2S)) 



= 2.5 ±0.6. (19) 



2.8. (20) 



Br(xw(2P)^ 7 T(lS)) 

One can see that it's agreeable with the experimental value. The full decay widths of Xhi(lP) and Xfci(2P) can 
be estimated by the branching ratios and our predicted decay widths. The results are: r XM (ip) <~ 85.7 keV, and 
F Xbl ( 2 P) ~ 66.5 keV. 



9 

In conclusion, we first calculated the radiative El decay width of Xci(lP)- The excellent agreement between our 
result and experimental value shows that this method we used is good to deal with the charmonium radiative decays. 
Then with the traditional radial excited charmonium state Xci(2P) assignment for X(3872) we calculated the radiative 
El decay widths of this particle, T(X(3872) -> 7J/VO = 33.0 kcV, T(X(3872) -> 7^(25)) = 146 kcV and T(X(3872) -> 
7^(3770)) = 7.09 keV. The value of r^/ 7 /r^, 7 is 4.4, which is consistent with experimental result by BaBar, but is 
larger than the up-bound reported by Belle recently. 

We also estimated the radiative El decay widths of the bottomonia states xm(1P) an d Xm(2P). Results are 
r(xbi(lP) -> 7T(1S)) = 30.0 keV, r(x 6 i(2P) -> 7T(1S)) = 5.65 keV, and r(x&i(2P) -> 7 T(2S)) = 15.8 keV. The 
predicted ratio Fy' 7 /Fy 7 of Xbi(2P) is consistent with experimental data. The full decay widths of r(x&i(lP)) = 85.7 
keV and r(x&i(2P)) = 66.5 keV (by the channel Xbi(2P) — > 7T(1S)) are also estimated. 
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